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Dynamical instabilities at the interface between two Bose-Einstein condensates that are moving 
relative to each other are investigated using mean- field and Bogoliubov analyses. Kelvin-Helmholtz 
instability is dominant when the interface thickness is much smaller than the wavelength of the 
unstable interface mode, whereas the counter-superflow instability becomes dominant in the opposite 
case. These instabilities emerge not only in an immiscible system but also in a miscible system 
where an interface is produced by external potential. Dynamics caused by these instabilities are 
numerically demonstrated in rotating trapped condensates. 



I. INTRODUCTION 

When wind blows over a water surface, the relative mo- 
tion between the water and the air generates instabilities 
at the interface; these instabilities in turn produce waves. 
Such interfacial instabilities between two fluids that are 
moving relative to each other are referred to as Kelvin- 
Helmholtz instabilities (KHIs) p]-[3| and they are found 
throughout nature. A system of superfluids is an ideal 
testing ground for KHIs because of the absence of vis- 
cosity. The Helsinki group experimentally realized shear 
flow between the A and B phases of superfluid 3 He in 
a rotating cryostat. They observed vortices penetrating 
from the A phase into the B phase due to KHIs 0, [H[ . 

Recently, KHIs and nonlinear dynamics of quan- 
tized vortices have been investigated for a gaseous two- 
component Bose-Einstein condensate (BEC) with rela- 
tive velocity [6]. In Ref. [fj, it is assumed that the two 
components are strongly segregated and that the inter- 
face thickness is much smaller than the wavelengths of 
unstable interface modes. If, however, the two compo- 
nents are weakly segregated, they strongly overlap with 
each other at the interface, at which the two components 
coexist with a relative velocity. Law et al [7[ showed that 
two miscible BECs moving through each other exhibit a 
dynamical instability, which we refer to as a counter- 
superflow instability (CSI). The dynamics generated by 
CSIs have recently been experimentally observed [8]. We 
therefore expect that KHIs compete with CSIs as the 
thickness of the interface is increased and that CSIs be- 
come important in the weakly segregated regime. 

In the present paper, we investigate both KHI and CSI 
on an equal footing in a phase-separated two-component 
BEC with relative velocity. We show that KHIs and CSIs 
are respectively generated for thin and thick interfaces; 
the interface thickness is controlled by the repulsive inter- 
action between the two components and external forces. 
The numerically obtained stability boundaries in the pa- 
rameter space are well described by the KHI in the thin 
interface limit and by the CSI in the uniform overlap 
limit. We propose various experimental situations for 



observing these instabilities in trapped systems. 

This paper is organized as follows. Section ITT1 reviews 
KHIs and CSIs in BECs, and performs the mean-field 
and Bogoliubov analyses for a system with a flat inter- 
face. Sections IIIII and IIVI propose experimental systems 
for generating these instabilities in pancake-shaped and 
cigar-shaped traps, respectively, and numerically demon- 
strate the nonlinear dynamics caused by these instabili- 
ties. Section [V] presents the conclusions of this study. 



II. INSTABILITIES AT AN IDEAL INTERFACE 

A. Kelvin-Helmholtz instability 

We first formulate a KHI in a strongly segregated 
two-component BEC with relative velocity between its 
two components [9[. In the mean- field theory, the two- 
component BEC is described by the macroscopic wave 
functions ipi and ip2- The Lagrangian for the system is 
given by 



L = 



J dr (Pi +P 2 -#12 |^i I 2 1^2 1 
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with rrij and Uj being the atomic mass and the exter- 
nal potential of the jth component, respectively. The 
inter- and intra-component interaction parameters have 
the form, 



gjj> = 2tt% 2 a jj'{m- 1 + m-, 1 



(3) 



where a^> is the s-wave scattering length between the 
jth and j'th components. The interaction parameters are 
assumed to satisfy the phase-separation condition [To| , 



#11#22 < 012- 



(4) 
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We assume that components 1 and 2 are respectively 
located in y < and y > and that the interface between 
the two components is located near the y = plane. In 
this subsection, we neglect the interface thickness and 
introduce an interface tension coefficient a [lll-ll3| , which 
originates from the excess energy at the interface. The 
Lagrangian can thus be rewritten as 

// pn poo \ 

dxdz dyPi + J dyP 2 j - aS, (5) 

where y = rj(x, z, t) is the position of the interface and 
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is the area of the interface. Taking the functional deriva- 
tive of the action J dtL with respect to r)(x, z, t) and set- 
ting it to zero, we obtain 



Pi(y = v) - ^2(y v) 



( d 2 r] d 2 r] 
V dx 2 dz 2 



0, 



(7) 



which corresponds to the Bernoulli equation in hydrody- 
namics. 

We consider a stationary state in which the jth com- 
ponent flows in the x direction with a velocity Vj as 



n o(y) ex P 



- (-fijt + rrijVjx) 



(8) 



where pj is the chemical potential for the jth component. 
The potential Uj(y) and the density distribution rij(y) 
are assumed to depend only on y, where ri\ = for y > r] 
and n 2 — for y < r] are satisfied. Substituting Eq. (|SJ) 
into Eq. (0) with r] = gives the equilibrium condition 
for the pressure, gnn 2 /2 = ^22^2/^- 

We assume that the system is approximately incom- 
pressible. A small deviation from the stationary state in 
this case is described as 
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vE^-exp iAje ( iyky cos(kx - uot) 
a sin(/cx — ojt), 



(9) 
(10) 



where Aj and a are infinitesimal parameters. From 
the kinematic boundary condition, the interface veloc- 
ity in the y direction (d/dt + Vjd/dx)r] must be equal to 

h/ (imjUj)^* dipj / dy\ y=r] , giving 

- (-l)J A jke-^^ = (Vjk-u)a. (11) 
rrij 

Substituting Eqs. ([9])- (fTT]) into Eq. ([7j) and neglecting 
second and higher orders of Aj and a, we obtain 



Pi 
k 



{w - V\k) 2 - fin. 



si 



P2 



(uj - V 2 k) 2 - f 2 n S 2 + ak 2 , 
(12) 



where n sl = m(rj + 0_), n s2 = n 2 (r] + 0+), pj = mjn sj , 
and fj = Uj(rj). Equation ([T2]) gives the dispersion rela- 
tion, 



(piVi + P2V2)k 
Pi + P2 



±1 



PiP2{V 1 -V 2 ) 2 k 2 , Fk 



yk 3 



(P1 + P2) 2 



pi 



- P2 

(13) 

where F = n s i/i - n s2 /2- 

We note that Eq. (|13p has the same form as the dis- 
persion relation for the KHI in classical incompressible 
and inviscid fluids [3]. When F > and 



(Vl - V 2 f > 4aF 



(Pi + P2? 
P\P\ 



V 4 



(14) 



a dynamical instability arises for fc_ < k < k+ : where 
PiP2(V 1 -V 2 ) 2 



k± 



1 ± 



V 4 

v cr 

(Vi - v 2 y 



2a(pi + P2) 
The most unstable wave number is given by 
, PiP2(V 1 -V 2 ) 2 



3a(p! + p 2 ) 



3Kt 



4(Vi - V 2 ) 4 



(15) 



(16) 



For |Vi — V2I < Kr, the system is dynamically stable. 
When F = 0, vanishes and the system is dynamically 
unstable for \V\ — V 2 \ > 0. In this case, the range of 
unstable wave numbers is given by 



< k < 



PlP2(Vl-V2) 2 

a(pi + p 2 ) 



(17) 



From Eq. (p~3]) . we find that the group velocity duo/dk 
diverges for k — >• and F > 0. This unphysical behavior 
originates from the incompressibility approximation in 
Eq. (|9]). For F < 0, the system is always dynamically 
unstable even when V\ — V 2 
a Rayleigh-Taylor instability jl4l4l6j. 



- 0. which is referred to as 



B. Counter-superflow instability 

We review the derivation of the Bogoliubov spectrum 
for a system of a uniform two-component BEC with a 
relative velocity [7|. The functional derivative of J dtL, 
where L is the Lagrangian in Eq. (pQ), with respect to ipj 
gives the Gross-Pitaevskii (GP) equation (j ^ j 7 ), 

dijj 1% 2 

(18) 

We assume Uj = and the miscible condition, 

#n#22 > g\ 2 , (19) 

in this subsection. We consider a small excitation <pj 
above a uniform state with a velocity Vj — hKj/rrij as 

i/jj = (v^ + <t>j) e- i ^ t ^ h+iK ^ r , (20) 
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where /ij = gjjfij + g 33 >ny + ft K 2 /(2rrij). Substituting We expand the small excitation as 
Eq. ([20]) into Eq. ([15]) and taking the first order of </>j, 
we obtain (j ^ j') 



zft- 



St 



ft 

(V + zK,) 2 - ^ + 2<? ii n i + fcvn^ 



2m 



ik-r — iojt 



v jk e 



-ik-r+icot 



(22) 



(21) and substitute it into Eq. ([2T]) . giving (j ^ j') 
I 



2m, 



(A: 2 + 2fc • Kj) + 



2m, 



ifc - QjjKjVjk + 933'\J n 3 n 3' ( u j'k Vj f k) = flUJU jkj 



(k 2 - 2k • 



v jk - gjjfijUjk - gjj'y/njrij' (u rk - v fk ) 



huj = ^(Vi + V 2 ) • k ± el + e 2 + 2e u ± 2 (ege? + 2e eN + ^i 



2„.2 \V2" 
2j 



1/2 



(23a) 
(23b) 



Diagonalizing the eigenvalue equation ([23]) . we obtain the Bogoliubov excitation spectrum. 

For simplicity, we assume that mi = rri2 = m and g\\n\ = #22^2 = u. Then, the eigenvalue of Eq. ([23]) has a simple 
form 



(24) 



where eq = ft 2 /c 2 /(2m), e r = Tik • (Vi — ^)/2, and U12 = #12 (^1^2) 1 ^ 2 - The expression in the square brackets in 
Eq. ([24]) for the negative sign becomes negative for 



max[fc r 2 cos 2 x ~ 4m(w + i/^/ft 2 , 0] < k 2 < max[/c 2 cos 2 x — 4m(w — ^i2)/ft 2 , 0], 

I 



(25) 



where the function max yields the maximum value of the 
arguments, x is the angle between Vi — V2 and fc, and 
fc r = m| Vi — V^l/ft. From Eq. ([25]) . the system is dynam- 
ically stable when A: 2 < 4m (u — ^12) /ft 2 . Minimizing the 
expression in the square bracket of Eq. ([24]) with respect 



to k and we find cos x — lj thus, the most unstable 
wave vector is parallel to the relative velocity. 



Crossover between Kelvin— Helmholtz and 
counter-superflow instabilities 



In Sec. Ill Bl we showed that a CSI emerges when the 
two components overlap with a relative velocity. We ex- 
pect that a similar situation arises when the interface 
between the two components is sufficiently thick in the 
system discussed in Sec. Ill Al since in the interface re- 
gion the two components overlap considerably and have 
a relative velocity. In this subsection, we show that a CSI 
emerges when the interface thickness is much larger than 
the characteristic wavelength. We also show that similar 
behavior is observed in miscible condensates separated 
by a potential gradient. 

We consider three types of interfaces: thin and thick 
interfaces produced by the intercomponent repulsion and 
an interface produced by a potential gradient; these in- 
terfaces are shown in the top panels of Figs. [U (a), [I] (c), 
and[T](d), respectively. The interface thickness produced 



by the intercomponent repulsion is characterized by 

9 



(26) 



where g = gu = #22 • For A -C 1, the density distribution 
is approximated as [12] 



n 



n i(y) = tt 



(-lytanh- 



/ 2&y 



(27) 



where £ = ft/ '(mgno) 1 / 2 is the healing length. The in- 
terface thickness is therefore £A -1 / 2 . The potential 
gradient in Fig. [T](d) has the form U\ — —U2 = Gy with 
G > being a constant, which stabilizes the interface at 

y = 0. 

To observe KHIs and CSIs in a phase-separated BEC, 
we perform Bogoliubov analysis. Substituting the sta- 
tionary state (jSJ) and a small excitation of the form, 



= u jk x (y) e 



ik x x — iujt 



3 K '- 



(y) e 



-ik x x-\-iujt 



(28) 



into the GP equation and taking the first order of <pj , we 
obtain the Bogoliubov-de Gennes (BdG) equation simi- 



lar to Eq. ([23]) . in which Kj has only the x component 
and rij, Ujk, and Vjk are replaced by rij(y), Ujk x (y), and 
Vjk x (y), respectively. The BdG equation is diagonalized 
numerically. 

The middle panels in Fig. [1] show the maximum value 
of the imaginary part of the Bogoliubov frequency Im uj 
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FIG. 1: (Color online) Density distributions of the ground states (top panels), maximum values of the imaginary part in 
the Bogoliubov spectra (middle panels), and profiles 5m(x,y) and 5ri2(x,y) of the most unstable modes defined in Eq. ([34)) 
(bottom panels). The field gradient G is in (a)-(c) and G — m^G/Ti 2 — 10 -7 in (d). In the middle panels, the solid lines 
indicate the boundaries between dynamically stable and unstable regions and the dashed lines are plots of Eq. ([29]) in (a) and 
Eq. ([31[) in (c). The color bars are normalized by gno. In the bottom panels, the relative velocity is V T = m^V Y /h = 0.08 and 
the field of view is 1600£ x 800£. All the results are independent of Vi + Vi- 



for each relative velocity V r = V\ — V2 and wave num- 
ber k x of the excitation, where the solid lines divide the 
dynamically stable and unstable regions. The stability 
boundary for the KHI in Eq. (fT7|) with p\ = p 2 = mno 
has the form 



mnoVj 2 
2a 



(29) 



The dashed line in the middle panel of Fig. Q] (a) plots 
Eq. ([29]) with the interface tension coefficient in Refs. [HI, 
Ell as 



Tin, 



3/2 



(30) 



The solid and dashed lines agree well for a low relative 
velocity. If G > 0, the stability boundary deviates from 
Eq. fl2SD for k H> @. The dashed line in the middle 
panel of Fig. [1] (c) plots 



m 



(31) 



which is the stability boundary for the CSI in Eq. ([25]) 
with u = u\2 and X = 0. The dashed line almost overlaps 



with the solid line. These facts indicate that the KHI is 
the dominant instability for the parameters in Fig.[T](a), 
whereas the CSI dominates in Fig. Q](c). 

The dominant instability changes from the KHI to the 
CSI when the interface thickness or the relative velocity 
increases. The crossover between the two instabilities can 
be estimated by equating Eqs. ([29]) and ([3T]h 



mnoV r 2 
2a 



which can be rewritten as 

h £ 

mVr 



(32) 



(33) 



using Eq. ([30]) . For example, for A = 2 x 10 3 , the 
crossover velocity is V r = m^V T /?t ~ 0.04. This is con- 
sistent with the middle panel of Fig. [T] (b), in which 
the stability boundary is parabolic for V r <C 0.04 and 
linear for V r ^> 0.04. Equation ([33]) indicates that the 
KHI is the dominant instability when the interface thick- 
ness ~ ^A -1 / 2 is much smaller than the characteristic 
wavelength h/(mV r ) associated with the relative veloc- 
ity. The CSI is dominant in the opposite limit, i.e., 



5 



£A -1 / 2 > h/(mV r ). It should be noted that the sta- 
bility boundaries in the middle panels of Figs. Q] (a)-Q] 
(c) have the universal form scaled by A 1 / 2 . This is be- 
cause the healing length £ is not relevant in the present 
problem, and the characteristic length scale for G = is 
only the interface thickness £A -1 / 2 . In fact, plotting the 
boundaries with respect to k x ^/S}^ 2 and V^A 1 / 2 , we find 
that they agree very well. 

Another difference between the KHI and the CSI is 
found in the excitation modes. The bottom panels in 
Fig. Q] show the change in the density by excitation of 
the most unstable mode, 



(a) A=10 



Snj(x,y) 



yjnj(y) + cu jkx (y)e l 



cv 



(34) 



where c <C 1 is a small constant. The first term on the 
right-hand side of Eq. ([34|) is the density distribution 
with a small excitation of the mode. We see that the 
excitation modes are localized near the interface. In the 
bottom panel of Fig. [T] (a), the wavelength of the mode 
is larger than the interface thickness, and the excitation 
of the mode shifts the interface sinuously. On the other 
hand, in the bottom panel of Fig. [T] (c), the wavelength 
of the mode is smaller than the interface thickness, and 
the dynamical instability develops stripes in the interface 
region. Interestingly, the stripes are inclined with respect 
to the interface, whereas the most unstable wave vector 
is parallel to the relative velocity in the uniform system 
discussed in Sec. Ill Bl 

Figure [T] (d) shows the case of the interface produced 
by the field gradient G > 0, where the miscible condition 
A < is satisfied. The stability boundary in the middle 
panel of Fig. [T] (d) is similar to that for the immiscible 
case with the field gradient (Fig. 2 of Ref. [6[). The exci- 
tation profile in the bottom panel of Fig. [T] (d) is inclined 
stripes, just as in Fig. [T] (c). If we apply the same field 
gradient to the system of Fig. Q] (c), we obtain results 
similar to Fig. [T] (d) (data not shown). 

Figure [2] demonstrates the time evolution of the im- 
miscible system for A = 10 ~ 2 and 10 _5 with no field 
gradient G = and the miscible system for A = — 10 ~ 5 
with a field gradient G > 0, where the relative veloc- 
ity is V r = 0.1 with V\ + V2 = 0. The initial state is 
the stationary state of the GP equation with small white 
noise to break the translation symmetry and trigger the 
dynamical instability. The time evolution is obtained by 
solving the GP equation by the pseudospectral method. 
For A = 10~ 2 [Fig. [2] (a)], the KHI is dominant, corre- 
sponding to Fig. [I] (a). We see that the interface is modu- 
lated due to the KHI (t = gn t/h = 5000). In Figs. El (b) 
and [2] (c), the interface is much thicker and the CSI is 
dominant; modulation arises in the region where the two 
components overlap. The nonlinear dynamics and vortex 
creation are quite different in the cases of the KHI and the 
CSI. In the former, quantized vortices are formed at the 
peaks and troughs of the wave, and enter into each com- 
ponent from the interface [6]. In contrast, in the latter, 
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FIG. 2: (Color online) Time evolution of the two- component 
BEC for (a) A 10" 2 with G = 0, (b) A = 10" 5 with (5 = 0, 
and (c) A = — 10 -5 with G = 10 -7 . The relative velocity is 
V Y = 0.1 with V1+V2 = 0. The four upper rows are the density 
profiles l^i] 2 and the bottom panels are the phase profiles arg 
ipi. The arrows in the bottom panel indicate an example of 
a vortex pair. The density is normalized by grtQ. Time t is 
normalized as t = gnot/h. The field of view is 1000£ x 750£. 



vortex pairs are created in the interface region (arrows 
in the bottom panels), which are subsequently disturbed 
in a complicated manner. In all the cases in Fig. [2j the 
total density \^i\ 2 + 1^1 2 is almost constant throughout 
the dynamics. 



III. INSTABILITIES IN PANCAKE-SHAPED 
SYSTEMS 

A. Kelvin-Helmholtz instability 

We have considered an ideal flat interface in Sec. |TT1 
In the following sections, we propose realistic systems 
for experimental observation of interfacial instabilities in 
trapped BECs. 

We first consider an axisymmetric harmonic potential 
given by 



U^r) = U 2 (r) = \m[Lol{x 2 + y 2 ) + w 



2_2i 



(35) 



where the radial and axial trap frequencies are uj± = 
2tt x 80 Hz and uj z = 2tt x 4 kHz, respectively. Since 
huj z is much larger than other characteristic energies, 
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FIG. 3: (Color online) Time evolution of the density and 
phase profiles. In the initial state, component 2 rotates with 
vorticity 712 = 18 and component 1 is at rest. The total num- 
ber of 23 Na atoms is N = 1.63 x 10 6 with an equal population 
in each component. The field of view is 93 x 93 /mi. The unit 



of density is N/a ± with a± = [h/(mu± 



il/2 



we reduce the system to two dimensions (2D) in the 
simulation. Assuming that the wave function can be 
written as ip(r) = ifj±(x^y)ifj z (z) with ip z (z) being the 
ground state of the harmonic oscillator, and integrat- 
ing the GP equation with respect to z, we obtain the 
2D GP equation with the effective interaction coeffi- 
cient gjjt = [muj z /(27rh)] 1 ^ 2 gjj' . Components 1 and 
2 are assumed to be respectively the hyperfine states 
\F,m F ) = |1,0) and |1, 1) of 23 Na atoms. The scattering 
lengths measured in Refs. [I?], G3 give an ~ 53.4<2b and 
&12 = «22 — 54.2<2b with ae being the Bohr radius, which 
satisfy the condition of phase separation (j4]). A strong 
magnetic field suppresses the spin exchange dynamics, 
1 0, 0) 1 0, 0) — >• |1,1)|1,— 1), due to the quadratic Zeeman 
effect jl9j . 

The initial state is prepared as follows. We first cal- 
culate the ground state with a centrifugal potential 
ft 2 n 2 I '[2m(x 2 + y 2 )] by the imaginary-time propagation 
of the GP equation. We then give the phases as 



tl>j 



Vil (7 = 1,2), 



(36) 



where the integer rij is the vorticity of the jth compo- 
nent and = arg(x + iy). This procedure is efficient 
for preparing an axisymmetric rotating stationary state. 
We thus obtain a stationary state with shear flow at the 
interface, as shown in the top panels in Fig. [3j where 
ni = and ri2 = 18. Because CL22 > Q>n and the cen- 
trifugal force on component 2, component 1 occupies the 
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FIG. 4: (Color online) Time evolution of the density profile 
of component 2 for 712 = 12. Other parameters are the same 
as those in Fig. [3] 



center and component 2 surrounds it. The presence of 
component 1 stabilizes the vortices in component 2. We 
add small numerical noise to the initial state to break the 
axisymmetry. The two components are assumed to have 
the same number of atoms in the following analysis. 

Figure [3] shows the subsequent time evolution. At 
t = 40 ms, a wavy pattern develops at the interface due 
to the KHI [6], which has approximately nine- fold sym- 
metry. The wave at the interface then grows (t — 44 ms) 
and quantized vortices are released from the peaks and 
the troughs of the wave into both components (t = 50 
ms). This interface behavior is similar to that in the flat 
interface [Fig. [2] (a)]. When the relative velocity at the 
interface is below a critical velocity (77,2 < 12), symmetry- 
breaking dynamics does not occur, whereas the flat in- 
terface is always unstable for a nonzero relative velocity 
[see Eq. (fT7|) ]. This indicates that in the trapped sys- 
tem the effective force is exerted on the interface [see 
Eq. ([14|)], which originates from the trap potential and 
the centrifugal force. 

Figure [H shows the time evolution of component 2 for a 
relative velocity just above the critical value (72,2 = 12). A 
five- fold pattern first emerges [Fig. |4] (b)], which changes 
to a four-fold pattern [Fig. |4] (d)] and then returns to 
a five-fold pattern [Fig. |4] (f)]. In contrast to the case 
of ri2 = 18 shown in Fig. [3j the deformation of the in- 
terface is moderate and no quantized vortices are gener- 
ated throughout the dynamics. Estimating the unstable 
wavelengths in Figs. [3] and |4] from the analytic expression 
(fTB|) is difficult because of an ambiguity in the interface 
tension coefficient a in inhomogeneous systems with a 
curved interface. 
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B. Bogoliubov analysis 



y 2 ) 1 / 2 . We expand the small excitation <j)j as 



We perform Bogoliubov analysis of a quasi-2D trapped 
system. We write the wave function as 

Mr, 0) = [fj(r) + 4>j{r,6)] e^e"^, (37) 



4>j(r, 9) = u jL {r)e iL0 -^ - v* L (r) e - iL9+iut , (38) 

where L is an integer and the excitation has L-fold sym- 
metry. Substituting Eqs. ([37|) and ([38]) into the GP equa- 



where fj(r)e inj is the stationary state, <pj is a small tion and taking the first order of 0j, we obtain the BdG 
excitation, \±j is the chemical potential, and r — (x 2 + equation (j ^ j'), 
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FIG. 5: (Color online) Imaginary part of Bogoliubov fre- 
quency uj = uj/uj± for m = and ri2 — 12 (circles), 
ri2 = 15 (squares), and ri2 = 18 (triangles). The excitation 
mode has L-fold symmetry. The total number of atoms is 
N — 1.63 x 10 6 . The conditions for triangles and circles are 
the same as those in Figs. [3] and [H respectively. The solid 
arrows indicate the modes relevant to the patterns in Figs. [3] 
and [H The inset shows a snapshot of the density profile of 
component 2 obtained by time evolution for 712 = 15. 



We numerically diagonalize Eq. (|39|) . 

Figure [5] shows the imaginary part of the excitation 
frequency for the same conditions as those in Fig. [3] 
(ri2 = 18, triangles) and Fig. |4] (77,2 = 12, circles). For 
ri2 = 18, the imaginary part is a maximum for L = 9, 
which is consistent with the nine-fold pattern in Fig. 03 
The emergence of the four and five- fold patterns in Fig. |4] 
is understood from the fact that the corresponding imag- 



inary parts are close to each other (the circles indicated 
by the arrows in Fig. [5]). In fact, we numerically con- 
firmed that for ri2 = 18, a nine-fold pattern or a 10-fold 
pattern emerges depending on the initial random noise. 
For ri2 < 12, no imaginary part appears. For the condi- 
tion in Fig.[5j only the lowest mode becomes dynamically 
unstable for each L. 



C. Counter-superflow instability 

We next consider the CSI in the pancake-shaped sys- 
tem. We use the hyperfine states |1, 1) and |1, 0) of 87 Rb 
atoms for components 1 and 2. According to Ref. [20|, 
the scattering lengths are an = a\2 = 100. 4ae and 
«22 = 100. 9<2b, which satisfy the miscible condition ([T9]) . 
We employ the harmonic potential plus a central optical 
plug given by 

1 



U 1 (r) = U 2 (r) = -m[u{(x 2 + y 2 )+uj 2 z z 2 



-f3(x 2 +y 2 )/< 



(40) 



where a and j3 are dimensionless parameters respectively 
characterizing the strength and the width of the optical 
plug, and a± = [h/ (moj^)] 1 / 2 . The same trap frequencies 
as those in Sec. lIII Al are used and hence the calculation is 
performed in 2D. We choose the parameters of the optical 
plug to be a = 1000 and j3 = 0.01. The optical plug is 
necessary for m ^ and ri2 ^ to eliminate undesirable 
instability. If the optical plug is absent, one component 
flows into the vortex core of the other component and 
the axisymmetry is broken before the CSI emerges. 

Figure [6] shows the time evolution of counterrotating 
condensates with initial vorticities n\ = —20 and 712 = 
30. Despite the miscible condition, the two condensates 
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FIG. 6: (Color online) Time evolution of the density profiles 
for counter-rotating miscible condensates. The vorticities of 
the initial state are m = —20 and ri2 = 30. The total number 
of 87 Rb atoms isiV = 3.19xl0 6 with an equal population in 
each component. The field of view is 84 x 84 fim. 



are weakly separated (t = ms of Fig. [6]) because of the 
difference in the centrifugal force. As time develops, a 
saw- toothed pattern emerges at the interface, as shown 
in Fig. [6l This behavior at the interface is similar to that 
in Figs. [2] (b) and [2] (c), indicating that the symmetry- 
breaking dynamics is driven by the CSI. The relation 
between the CSI and quantum turbulence has recently 
been studied [211 ]. 



IV. INSTABILITIES IN CIGAR-SHAPED 
SYSTEMS 

A. Kelvin— Helmholtz instability 

We propose two systems for observing the KHI in a 
cigar-shaped geometry. 

The first one is similar to that in Sec. IIII A( non- 
rotating component 1 is surrounded by rotating compo- 
nent 2. We use a harmonic potential in Eq. ([35]) with 
ujj_ = 2ir x 40 Hz and uj z = 2tt x 8 Hz. Components 
1 and 2 are respectively the hyperfine states |2, 1) and 
1 1,-1) of 87 Rb [22[. The ratio between the scattering 
lengths of these states is an : an : a22 = 0.97 : 1.0 : 1.03 
with their average being 5.5 nm [23|, [24|, which satisfy 
the immiscible condition (jl]). 

We solve the GP equation by performing a 3D numer- 
ical simulation. Figure [7] (a) shows the time evolution of 
the cross-section densities \^\{x, y, 0)| 2 and |t/?2 (#,?/, 0)| 2 
for the initial vorticities rt\ = and ri2 = 7. The ini- 
tial state has an axisymmetric interface with a relative 
velocity. At t = 72 ms, a wavy pattern with four- fold 
symmetry develops at the interface due to the KHI. Vor- 
tex lines are then released from the interface into both 
components at t = 84 ms as shown in Figs. [7] (a) and [71 
(b). For a smaller relative velocity at the interface (e.g., 
774 = and 77,2 = 1), symmetry breaking dynamics does 
not occur, indicating the existence of a critical velocity, 



(a) component 1 component 2 (b) 



ms 



72 ms 



84 ms 





o 


O 


o 


tft 


© 


■ 





FIG. 7: (Color online) (a) Time evolution of the cross-section 
densities y, 0)| 2 (left panels) and \ip2(x, y, 0)| 2 (right 

panels) for m = and 772 = 7. The field of view is 20 x 20 
/mi. The total number of 87 Rb atoms is N = 1.75 x 10 5 with 
an equal population in each component. The bottom panels 
show the phase profiles at t = 84 ms. (b) Isodensity sur- 
faces of component 1 (inside) and component 2 (outside) at 
t = 84 ms. For clarity, only the surface of component 2 is 
transparent. 



as in the case of Sec. IIII A[ 

The second configuration for observing the KHI in a 
cigar-shaped trap is shown in Fig. [8] (a), where compo- 
nents 1 and 2 are located in the z < and z > regions 
and their interface lies on z ~ 0. We use the same hy- 
perfine states of 23 Na as those used in Sec. IIII Al The 
potential in Eq. (|ID|) is used with uj± = 2tt x 40 Hz, 
uo z = 2tt x 8 Hz, a = 50, and /3 = 0.4. The optical plug is 
applied to prevent one component flowing into the vortex 
core of the other component. We apply a magnetic field 
gradient of dB/dz = —17.3 mG/cm, which pushes com- 
ponent 2 (|1, 1)) in the +z direction. This force controls 
the strength of the KHI. 

Figure [8] (b) depicts the time evolution of the system 
with initial vorticities of 77 1 = 2 and 772 = —2. The 
relative velocity at the interface induces the KHI, which 
breaks the axisymmetry and deforms the interface [280 
ms of Fig. [8] (b)]. A quantized vortex is then released 
from the interface into each component at t ~ 298 ms. 
The vortex lines lie in the radial direction; in contrast, 
the vortex lines lie along the z direction in Fig. [TJ 



B. Counter-superflow instability 

As shown in Sec. [Ill the inter facial instability is dom- 
inated by the CSI when the overlap between the two 
components is large. We study this situation for a cigar- 
shaped system. The two components are the same hy- 
perfine states of 87 Rb as those used in Sec. IIII C( and the 
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FIG. 8: (Color online) (a) Column densities f dz\ipj\ 2 (left 
panels, 40 x 40 jim) and f dy\ipj\ 2 (right panels, 40 x 132 /mi) 
of the initial state for n\ — 2 and ri2 — —2. The white arrow 
indicates the direction of the Stern-Gerlach force produced 
by the field gradient of dB/dz — —17.3 mG/cm. (b) Time 
evolution of the isodensity surfaces of components 1 (lower) 
and 2 (upper). The solid arrows at t = ms indicate the 
directions of the rotation. The phase profile of component 2 at 
t = 298 ms is shown in the rightmost image. The total number 
of 23 Na atoms is N = 6.6 x 10 5 with an equal population in 
each component. 

potential is the same as that in Sec. IIV Al 

Figure [9] shows the time evolution of the system for 
initial vorticities of rt\ = 3 and ri2 = —3. Since the 
scattering lengths satisfy the miscible condition and the 
field gradient is small, the two components in the initial 
stationary state widely overlap with each other [Fig. [9] 
(a)]. The time evolution in Fig. [9] (b) is similar to those 
in Figs. 121(b) and[2](c) in the sense that stripes of the two 
components are formed where they overlap. The spiral 
pattern in Fig. [9] (b) (t ~ 40 ms) corresponds to the 
inclination of the stripes in Figs. Q] (d) and [2] (c). 

V. CONCLUSIONS 

We have studied the dynamical instabilities at an in- 
terface of a two-component BEC with a relative veloc- 
ity between the two components. When the two compo- 
nents are strongly segregated and the interface thickness 
is negligible, a KHI is generated at the interface. On the 
other hand, when the interface thickness is much larger 
than the unstable wavelength, a CSI dominates over the 
KHI. We have proposed realistic experimental systems 




ms 40 ms 44 ms 



FIG. 9: (Color online) (a) Column densities f dz\ipj\ 2 (left 
panels, 27 x 27 /mi) and J dy\ipj\ 2 (right panels, 27 x 118 /mi) 
of the initial state for m = 3 and n<i — —3. The white arrow 
shows the direction of the Stern-Gerlach force produced by a 
field gradient of dB/dz = 0.34 mG/cm. (b) Time evolution of 
the isodensity surface of component 1. At t = 0, component 
2 is also shown as the transparent surface. The total number 
of 87 Rb atoms is N = 1.3 x 10 6 with an equal population in 
each component. 

of 23 Na and 87 Rb BECs in pancake-shaped and cigar- 
shaped traps. In a pancake-shaped trap, a KHI can be 
observed at the interface of two components separated 
inside and outside, and rotating with different vortici- 
ties (Figs. EH and H). The patterns produced by the KHI 
are consistent with the Bogoliubov analysis (Fig. EJ). Us- 
ing miscible two components, a CSI can also be realized 
(Fig. [6]). For a cigar-shaped trap, we proposed two con- 
figurations for observing the KHI: the two components 
separate radially (Fig. [7)) or in the z direction (Fig. [8]). 
The vortex lines released from the interface lie along the 
radial direction in the former, whereas they lie along the 
z direction in the latter. The interface thickness can be 
controlled using the miscible condensates with a mag- 
netic field gradient in the z direction (Fig. [9]). 

The initial states used in Sees. IIIII and IIVI are mul- 
tiply quantized vortices with different vorticities in the 
two components. Such states can be prepared using 
Laguerre-Gaussian beams [25| or topological phase im- 
printing [26] . Multiply quantized vortex states with large 
vorticities can in principle be generated by repeated ap- 
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